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Abstract 

Structural relaxation in hard-sphere colloidal glasses has been studied using confocal microscopy. 
The motion of individual particles is followed over long time scales to detect the rearranging regions 
in the system. We have used normal mode analysis to understand the origin of the rearranging 
regions. The low frequency modes, obtained over short time scales, show strong spatial correlation 
with the rearrangements that happen on long time scales. 
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Supercooled molecular liquids, polymers or metallic systems, when cooled or compressed 
in a way so as to prevent crystallization, undergo a glass transition. The transition to the 
glassy state is characterized by the slowing down of the dynamics of the system accompanied 
by a spectacular increase in the viscosity, to values on the order of ~ 10 13 Pa s pp. Both 
features are usually attributed to the fact that the motion of individual molecules or particles 
become arrested due to the presence of the cage formed by neighboring particles. However, 
on longer timescales the particles escape from their cages, leading to the relaxation of system. 

In recent experiments on colloidal glasses, the motion of individual particles was followed 
using confocal microscopy to investigate such cage rearrangement events in quiescent [21 E] as 
well as sheared [HE] systems. It was found that the particle dynamics is very heterogeneous 
in the sense that some regions exhibited much stronger activity than others: some small 
parts of the system are more susceptible to rearrangements than others. Earlier studies 
have reported such heterogeneous dynamics in glass forming liquids and polymers [6]. 

Perhaps one of the main questions about glassy dynamics is then how to understand 
and eventually predict the relaxation events that involve rearrangement of molecules or 
particles. One clue comes from recent computer simulation studies of supercooled liquids, 
which suggest that more susceptible regions result from localized 'soft' (i.e., low frequency) 
modes of the system [7]. The soft modes in glassy systems is a subject of much interest 
lately [8]- [15], as they are related to the anomalous low temperature properties of glasses. 
Using normal mode analysis, their existence has been reported recently in a few experiments 
on colloidal glasses [H] -[EE]- However, there is very little experimental evidence that connect 
these modes with the relaxation events in colloidal glasses |12j . 

In this Letter, we identify the low frequency modes in a hard-sphere colloidal glass while 
simultaneously studying the rearrangements. This is done by following the motion of flu- 
orescent colloids for a long time using rapid confocal microscopy. On the one hand, the 
particle motions can be analyzed in terms of the normal modes of the system. On the other 
hand, rearrangements can be identified by looking at the changes in neighbors. The main 
conclusion is that, indeed, the rearrangements happen along the softest modes of the system. 
We thus establish the structural origin of the rearranging regions by comparing them with 
the spatial maps of the amplitudes of the soft modes. 

We study the motion of particles in a dense hard sphere colloidal glass using confocal mi- 
croscopy. It is prepared by suspending sterically stabilized fluorescent polymethylmethacry- 
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late (PMMA) particles in a density and refractive index matching mixture of Cycloheptyl 
Bromide and Cis-Decalin. The particles have a diameter of a = 1.3/im, with a polydis- 
persity of 5% to prevent crystallization. The organic salt TBAB (tetrabutylammoniumbro- 
mide) is used to screen any possible residual charges. A dense suspension of volume fraction 
<p ~ 0.59 was prepared by diluting a sediment that was centrifuged to random close packing 
(^rcp ~ 0.64). Since the particles are poly- disperse, the volume fraction of the sample is 
subject to an error of ~ 2% [13]. 

Two dimensional slices through the 3d glass of fluorescent particles are acquired at a 
scanning speed of 10 frames per second using a fast confocal microscope (Zeiss LSM 5 live) 
in a field of view of lOO/xm x 100/im. The image plane was chosen so that it is sufficiently 
~ 25 away from the coverslip to avoid any possible boundary induced effects. Images 
are acquired for a period of total 1200s which results in 12000 independent snapshots. 
The positions of the particles are then identified in all the frames using standard particle 
tracking softwares and linked to construct the two dimensional trajectories. The mean square 
displacement (MSD) < 5r 2 > of the particles computed using these trajectories is shown 
in Fig. [TJ On short time scales, the rise of the average displacements is due to diffusion of 
particles before they start feeling the presence of their neighbors. Beyond this time scale, 
we observe a plateau in the MSD where the motion of a particle is restricted by a shell of 
nearest neighbors constituting the 'cage'. It extends roughly up to ~ 500s in the present 
experiments. Note that this times scale compares well with the previous measurements 
done on hard sphere colloidal glasses |3j. On even longer timescales, at the end of the 
plateau, the mean square displacement increases again. This corresponds to a long-time 
diffusive behavior known as the a-relaxation [161 E] regime and is usually attributed to 
cage rearrangements. The average displacement in the long-time diffusive regime is linear 
in time: < 5r 2 >~ t, as follows from Fig. [TJ the data in this region is indeed parallel to a 
line of slope unity. The effect of aging over the time scale of our measurement was found to 
be negligible. 

However what is not clear from these averaged quantities is the possible heterogeneity 
of the dynamics, as suggested by previous experiments [2] . The inset to Fig. [T] shows a 
single particle trajectory; there is typical 'cage rattling motion' followed for this particle 
by a rearrangement, in which the average of the center of mass changes positions. This 
does not happen for all particles over the observation time, and such events can be used 
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to see where rearrangements happen. To investigate this, we compute the square of the 
relative displacements df(Ai) of each particle with respect to its nearest neighbors over a 
time interval At from the reconstructed trajectories Xi(t),yi(t). The relative displacement 
df is defined as : 



where n denotes the number of nearest neighbors that is determined using the method of 
Delaunay triangulation and Ar^t) = — rj denotes the difference vector of the particles 
i,j at time t; this expression of df is in fact motivated by the non-affine measure of plasticity 
used in sheared amorphous solids [T§] . 

Figure |2^a) shows the probability distribution function of the relative displacements of the 
particles evaluated over a time scale of At ~ 1200s. The extended tail of the distribution 
implies the presence of large displacements in the system. To detect the highly mobile 
particles, we tag all those particles that undergo a relative displacement beyond a cutoff 
value r\. This cutoff r c = 0.245/im is ~ 1.25 times the square-root of the mean of the above 
distribution. The spatial distribution of the particles that have d\ > r\ is shown in Fig. [21(b) . 




Here, the red circles refers to the particles with large d\ (> r%) in a background of less active 
particles shown in pink. The particle with large df appear in the form of clusters in the field 
of view. Our observation shows that particles in these regions undergo significant change in 
their neighborhood. 

To demonstrate this, we identify the particles that change neighbors. The neighbors 
of each particle at any instant of time t can be identified using the method of Delaunay 
triangulation. We then compare the neighbors of each particle at time t and t + At. This 
gives us directly the number of neighbors lost or changed around a particle over a time 
interval At. In simulations of supercooled liquids [7], the irreversible rearrangements were 
identified as particles that have lost four neighbors. However, in experiments there are very 
few particles that lose four neighbors over the time scale of our observation. So, we identify 
those particles that have lost three or more neighbors in our data. Figure shows a plot of 
these particles (squares) along with the clusters of highly mobile particles. We observe that 
the majority of the particles that have lost three or more neighbors belong to the regions 
of large relative displacements. This establishes that df is a good way to detect rearranging 
regions in the system. 




(1) 
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Now that we have identified the 'active' regions, we will see whether the visualization of 
the low frequency modes allow us to predict where in the system rearrangements will take 
place. The existence of these modes was established by doing a normal mode analysis of 
particle trajectories over the plateau region of MSD [9]-[II]. There are hardly any rearrange- 
ments on this time scale, so we can reasonably assume that the particles are close to their 
local equilibrium, allowing us to determine the normal modes of the system. The displace- 
ment of the i th particle from the mean position is given by m = x~ < Xi >,Vi— < jji >. 
We then construct the covariance matrix of displacements [HI E] 

Dij = (uiUj), (2) 

where the matrix indices i,j = 1, 2, ...2N (N particles) runs over both the particle number as 
well as cartesian components in two dimension. The average () is done over a time interval 
At. By diagonalizing we obtain 2N eigenvalues Aj = 1, 2, 3..2N and their corresponding 
2N eigenmodes v(uii)s. The frequencies are related to the eigenvalues as u>i{= y/l/Xi). Note 
that the above method does not take into account any anharmonicity or damping effects that 
is present in the system [9]- [IT], [15]. The left panel of Fig. [3] shows the displacement vectors 
of one of the lowest frequency modes in the system that is obtained by averaging over a 
time interval At = 500s. The quasi- localized nature of the mode is apparent; some smaller 
regions have larger amplitude than the rest. The degree of localization of these modes can be 
measured by the average participation ratio following P{oj) = [N(J2i ^(^-^iO^)) 2 ] \ where 
Vi(u) is the normal mode (of frequency ui) amplitude projected onto the i th particle. For 
a strongly localized mode, where the whole amplitude is concentrated over a few particles 
in the system, P{ui) scales as 1/N, whereas for an extended mode it is of order of unity. 
The average participation ratio of the lowest frequency modes in our system that display 
quasi-localized character [3, [18] lies below P{ui) < 0.36. 

We will now use the spatial maps of particle participation ratio p%{ui) = vf |7J to visualize 
the localized nature of these modes. A contour plot of the particle participation ratios in 
Fig. [3] (right) illustrates the spatial structure of the mode shown on the left panel of the 
same figure. We can use these normal modes to predict the relaxation events provided 
that their structure remain invariant with time. To test the time invariance, we determine 
the low frequency modes by averaging over different intervals of time : 300s, 400s and 
500s. Figure [4] shows the contour plots of the particle participation ratios, averaged over the 
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25 lowest frequency modes, for different time intervals. We do not observe any significant 
change in the mode structures, besides slight variations in a few of the quasi-localized zones. 
This allows us not only to understand the heterogeneous dynamics but also to predict -based 
on short time measurements of the normal modes- where on long times the rearrangements 
are likely to happen. 

We therefore compare the contour map of particle participation ratio of the low frequency 
modes and the rearranging regions to establish a connection between the two. Figure [5] shows 
a superposition of the contour maps of the particle participation ratio, averaged over 25 
lowest frequency modes, along with the clusters of particles with higher mobility (Fig. [2]) in 
the field of view of 100 x lOO/xm 2 . Clearly, the 'active' regions where rearrangements happen, 
exhibit strong spatial correlation with the quasi-localized zones of the low-frequency modes; 
the rearrangements apparently originate in region of high particle participation ratio. It is 
worth pointing out that the normal modes were computed by averaging the particle motion 
over a time interval At = 500s, and the rearranging regions were detected in the later stages, 
over a time interval At = 1200s. The spatial distribution of the low frequency modes, 
thus represents the structure of the system much before the rearrangements occur. The 
robustness of our results were verified by performing several measurements. For example, 
see the supplementary figure that shows the contour map of particle participation ratio of 
the low frequency modes and the rearranging regions at <ft = 0.57. 

In conclusion, we have studied the structural relaxation of a hard-sphere colloidal glass. 
The dynamics of the particles exhibit spatial heterogeneity : clusters of highly mobile par- 
ticles are observed. The relaxation of the system occurs through rearrangement of particles 
in these highly active regions. To understand the nature and origin of these rearranging 
regions, we have computed the low-frequency modes of the system using normal mode anal- 
ysis. The low frequency modes display quasi-localized character, which is evident from the 
contour maps of the particle participation ratios. Our analysis reveals that the regions where 
rearrangements occur are spatially correlated with the quasi-localized zones of the low fre- 
quency modes. This demonstrates that the system indeed relaxes along the softest available 
modes. 

Even though these results are obtained for a colloidal glass under quiescent condition, 
they should be relevant to weakly driven colloidal glasses [3]. It is known for a long time 
that the plastic deformation in glasses occur in localized regions that are referred as shear 
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transformation zones (STZ) [U[19j|20]. However, their origin remains unclear. Our results 
should motivate further research along the direction of normal modes to identify the origin 
of the shear transformation zones. 
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FIG. 1: Mean square displacement (MSD) of the particles is shown as a function of time over the 
whole period of measurement ~ 1200s. The dark line has a slope of unity. Inset shows a typical 
trajectory a particle that undergoes cage rearrangement where average position shifts from one 
location to another. 




FIG. 2: Left: The probability distribution function of the relative displacements di for the present 
suspension measured over a time interval At = 1200s. The mean of this distribution is < df >= 
0.039/um 2 , and it indicated by the dashed vertical line. Center: A spatial map of the particles where 
the highly mobile particle with ad, 2 > r 2 (0.06/im. 2 ) are shown by the red circles along with the less 
active particles denoted by pink circles. Right: Superposition of the particles that have relative 
displacement df > r\ (circles) and the particles that have lost three or more neighbors (squares). 
Almost all the particles losing neighbors belong to the regions of larger relative displacement. 
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FIG. 3: Visualization of low frequency modes that are obtained by averaging the over a time 
interval At = 500s. Left: Example of one of the lowest frequency modes for the present system 
over the whole field of view 100 x 100/rni 2 . There are regions where the colloidal particles show 
higher activity than the rest in the field of view. The average participation ratio P(oj) = 0.17. 
Right: Contour plot of the local participation ratios for the same mode. 




At = 300 s At = 400 s At = 500 s 



FIG. 4: Contour plots of particle participation ratios, averaged over the lowest 25 normal modes 
in a field of view of 60/im x 60/im. This is shown for successively three different time intervals of 
averaging of Dij : 300s, 400s, 500s, from left to right. The mode structure remain largely invariant, 
besides slight variations in few of the quasi- localized zones. 
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FIG. 5: Superposition of the particle participation maps, averaged over lowest 25 modes, and the 
particles that have relative displacements beyond the cutoff r^, where r c = 0.245/im. The normal 
modes were obtained by averaging the Dij over a time interval At = 500s, and the rearranging 
regions were detected over a time interval At = 1200s. The rearranging regions display significant 
spatial correlation with the quasi-localized zones of the normal modes. 
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